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ABSTRACT 


Mimomeaiytreal anvestigation of laminar free convection 
from a semi-infinite vertical plate is undertaken with the 
following special conditions imposed. The plate is main- 
tained at a constant temperature greater than that of the 
ambient fluid. At an arbitrary distance above the leading 
edge a constant rate of blowing or suction is initiated 
extending uniformly along the remainder of the plate. The 
solution was obtained using the method of matched asymptotic 
expansions whereby an inner series described the velocity 
and temperature profiles near the discontinuity in the 
Vicinity of surface and an outer series approximated the 
behavior at greater distances from the plate. The inner 
series was also used to determine the heat transfer charac- 
teristics at the surface of the plate. The ambient fluid 


roseassumed to be air with a Prandtl number of 0.72. 
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I. INTRODUCTION 


ime study of free convection has been the subject of 
ongoing research over the past century both from an analyti- 
@al and experimental point of view. An excellent overview 
Seecme progress made in this area can be found in Ede [1]. 
Two outstanding expositions of the fundamentals of heat 
Meanster by the mechanism of free convection are presented 
by Gebhart [2] and Moore [3]. 

ie technigue of solving boundary layer problems by 
employing similarity analysis is discussed in general by 
Hansen |4| and specific applications of the method are pre- 
wemeead by Yang [5]. 

ime problem discussed in this paper is that of an analyti- 
meeminvestigation of the effects of blowing and suction 
applied at a specified distance above the leading edge of a 
Uniformly heated semi-infinite vertical plate on the heat 
transfer characteristics at the surface of the plate. 

An early study of the effects of mass transfer on free 
convection was conducted by Eichhorn [6] who showed that 
Bamtlar solutions are possible for blowing or suction rate 
distributions which vary as the distance from the leading 
Semenor the plate raised to a power related to the variation 


in surface temperature. 


1S) 





Pace the prescribed surface temperature itself varied as 
Mepewer Of the distance from the leading edge, this proved to 
be a somewhat restricted case of limited physical application. 
Mrenwcase Of a uniformly heated plate with constant blowing or 
mmerton alone the entire extent of the surface was succinctly 
treated by Sparrow and Cess [7] who extended their analysis 
to include independent power law variations of both surface 
temperature and mass transfer. Their analysis which was based 
Seam asymptotic series expansion of the stream function and 
temperature distribution evaluated the first order effects of 
blowing and suction on free convection. An extension of this 
work to include possible variations in the fluid properties 
was conducted by Parikh, Kays and Bershader [8] who employed 
meaumerical finite difference scheme to solve the resulting 
feeeerons, their numerical predictions were in close agree- 
ment with interferometric data obtained through their own 
experiments. The results of these authors and especially 
those reported by Sparrow and Cess have been of great value 
in affording a comparison with results obtained in the 
present study. A comprehensive investigation of the effects 
Seeesuction and blowing on heat transfer and wall skin fric- 
tion which included regions on the plate at large distances 
from the leading edge was conducted by Merkin [9]. In that 
analysis a step-by-step numerical solution of the boundary 
layer equations begins at the leading edge and proceeds 


Systematically up the plate until the asymptotic velocity 
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and temeperature profiles are obtained to whatever accuracy 
meerequired. 

There have been several additional studies in the area 
of free convection past a vertical plate with various condi- 
tions of wall temperature or heat flux and imposed variations 
mimeeranispiration at the surface of the plate. An analysis 
which transformed the laminar boundary layer equations re- 
mmecane from arbitrary distributions in both wall temperature 
and blowing into ordinary differential equations with variable 
coefficients and constant boundary conditions was presented 
by Vedhanayagan, Altenkirch, and Eichhorn [10]. In the area 
Of unsteady free convection where either the plate surface 
temperature or transpiration rate varies with time there have 
Been at least four notable studies during the past twenty 
years. Nanda and Sharma [11] undertook analytical inves- 
tigations and determined that similarity solutions were 
possible for the following two cases: surface temperature 
weoyame as some power of time and suction velocity varying 
as time to the minus one-half power; and the second case 
where surface temperature is an exponential or periodic 
function of time and the suction velocity is constant. The 
effects of mass diffusion on the unsteady free convective 
flow past a simi-infinite porous plate with constant suction 
were studied through mathematical analysis by Soundalgekar 
and Wavre [12]. With the assumption that the plate temperature 


oscillated in time about a constant mean, approximate 


eS 





solutions were derived for the mean flow, transient flow, the 
amplitude and phase of the skin friction and the rate of heat 
transfer. The effects of the nondemensional parameters Grashof 
number, Prandtl number, Eckert number and Schmidt number were 
also reported. The most general unsteady situation where both 
Pemeecemperature and suction velocity are periodic functions 
Beeecime was considered by Lal [13]. All of the previous 
memanmes Cited have neglected the effects of viscous dissipation; 
Memever, this effect is not insignificant when the natural 
memuectiOn flow field is very large or the flow occurs in the 
Deesence of extremely low temperatures or very high gravity 
meas. the combined effect of viscous dissipation with the 
feeeron velocity and plate temperature oscillating at the 
Seer requency on free convection past a vertical plate was 
Beoeeyezed by Soundalgekar and Pop [14]. In their study expres- 
SwomemtOr the mean velocity, transient velocity, transient 
Semeeratire profiles, amplitude and phase of skin friction 
Guemeene rate of heat transfer were derived and the results 
momeered tO those reported for the case of constant suction 
melocity. 

Peeeevirew Of the literature in the field of free convection 
oumeaevertical plate clearly indicates a desire on the part 
of researchers to investigate perturbations of the thermal 
and momentum boundary layers resulting from the introduction 
Seecemperature or velocity discontinuities along the surface 


seepene Dlate and their ultimate effect on skin friction and 
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Meee transier at the wall. The specific nature of the dis- 
continuity will determine the manner in which the perturbation 
peeeectS propagate through the boundary-layers and; therefore, 
will often indicate a preferred method of analytical solution. 
Although there are several possible mathematical techniques 
available to researchers in this field, two fundamentally 
Mieererent approaches predeminate, each with numerous variations. 
An example of the finite difference method of solution to the 
boundary-layer equations resulting from a discontinuity in the 
Sunbtace wall temperature can be found in Na [15]and Meena and 
Nath [16]. This technique has gained in acceptance with the 
meeert OL larger and faster digital computers. The other 
major technique involves the expansion of both the stream 
minetion and temperature profile in an asymptotic series. 
i@emonriginal proponent of this method was Goldstein who pre- 
Sented the details of the asymptotic series technique, 
commonly known as the Goldstein series, in [17]. After re- 
Meine the set of governing boundary-layer equations by the 


mmemeaniction of a stream function such that: 


BU ap 
oe 3y eiiiden Vie anaes 


two new variables are introduced which are of the following 
Pemeral form: 


1 
1) In the streamwise direction: €& = x? 


a2) 

n 

2) In the direction normal to the plate: n = yx 

A series is then formed in powers of & for the stream function 


and temperature distribution with undetermined coefficient 


stg 





functions fi respectively which depend only on 7. The 


2, 
meioem is then finally reduced to one of determining the 
Meerricient functions; the solution of a fifth order set of 
meear coupled ordinary differential equations for each pair 
@eeeoetticients subject to the appropriate boundary conditions. 
This technique is especially powerful when employed to inves- 
mmeaee CONditions at or near the surface of the plate. In 
order to examine the laminar free convection wake above a 
heated vertical plate Yang [18] successfully employed this 
technique to obtain detailed velocity and temperature profiles 
mimeene immediate vicinity of the trailing edge and then 
expanded these results to the rest of the wake by an integral 
Pemmtion. Another example of this technique is found in 
Kelleher [19] where a step discontinuity in the surface 
temperature imposed at a finite distance above the leading 
edge of the plate caused a perturbation of both the thermal 
and momentum boundary-layers. An inner expansion of the 
Goldstein type was employed to obtain the behavior of heat 
Peaster near the surface of the plate and an outer asymptotic 
series expansion was constructed to determine the velocity 
maceetemperature profiles at greater distances from the plate. 
A similar approach is employed to analytically investigate 

the present problem where a discontinuity in the velocity 
boundary condition is propagated through the momentum and 
thermal boundary-layers. Laminar flow is assumed as are 


eemsctant transport properties of the ambient fluid. The 
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brimary aim of the study is the determination of any en- 
hancement or degradation of the heat transfer rate at the 
surface of the plate compared to the same plate without 


Mowing Or SUCtion over the entire length. 
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Il. ANALYSIS 


See COORDINATE SYSTEM 

ine problem considered in this paper is that of laminar 
meee convective flow from a uniformly heated semi-infinite 
vertical plate with leading edge oriented downward in a 
Seanaard Gravitational field. A constant rate of blowing or 
Suction is applied at the surface beginning at an arbitrary 
o@ietcance, L, above the leading edge. A rectangular cartesian 
coordinate system with origin as indicated in Figure 1 was 
employed in the analysis. As is indicated throughout the 
Papen Various transformations of the basic coordinates are 


invoked to facilitate the solution of the problem. 


Pee GOVERNING EQUATIONS 

The analysis is based on the standard Boussinesq free 
convection approximation which assumes that only density va- 
Tiations giving rise to buoyant forces are considered with 
mitemertrtects Of viscous dissipation neglected. With these 
Standard assumptions the two dimensional equations of con- 


miniuicy, momentum, and energy boundary-layer form are: 


es 


Q? & 
alls] 
Qs 
< 
S 





Poet ad) (2) 
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Faoume io basac Coordinate System 


oo 








= T os 7 e 
u ua ey ms oe : (3) 
ax ay ay? 
where we 1s the constant ambient fluid temperature. These 


equations must also satisfy the following physical boundary 


Sonditions: 
Pitey = 0: u = 0 
WS es X>Q (32) 
=e Se) 
T = T 
W 
Moceas y + >: u = 0 
T=T 


20 
The non-zero normal velocity boundary condition at x = 0 and 
y = 0 introduces a discontinuity into the boundary-layer equa- 
tions which has a primary direct effect on the formation of 
the momentum boundary-layer and a secondary influence on the 
development of the thermal boundary-layer. <A positive value 
for = will represent the blowing of fluid away from the wall 
and a negative value Ore Wits Geprescnus fic sllellOn of .<t lurid 
through the plate. 

In order to facilitate the analysis, the variables in 
eam@ations (1)-(3) are non-dimensionalized according to the 


following transformations: 


(4a) 


c| 2s | 


— 
x = 
a 


(4b) 


mrs | 


eh 





— ul 
us = (4c) 
Ee 
7 Tae (4d) 
T- T, 
— (4e) 
=e 
W oO 


where L is the distance above the leading edge where suc- 
tion or blowing commences, and ne is the constant wall 

temperature. The equations are further simplified through 
the introduction of the conventional stream function w and 


Myestretching y, u, and v as follows: 


eos = y (2p * (5a) 
fe n(4Gr)* = 3 (5b) 
ve = v(64Gr)7* = - a [oe] 


where the Grashof number, Gr, is evaluated with respect to 
meeeene Length scale associated with the discontinuity. It 
should be noted that this length scale is not the conven- 
tional one used in free convection problems where distance 
poom tne Leading edge is normally used to evaluate Gr. Now 
meaeeresubsStituting equations (5) into equations (1)-(3), the 
Momtimuity equation is eliminated and the momentum and energy 


pauacions reduce to the following: 





+ aye ° oxdy® ~ + Ox ° By*? = yes * P (6) 





Q 
NM 
cD 


dW 90 Cas) 


ow : 
+ 3% Ox ~ + 3x By* (7) 


| 
ss 
Qo 
. 
ho 


The non-zero normal velocity component imposed at the surface 
of the plate perturbs the governing equations (6) and (7) 
Mmeeeluding the existence of a similarity solution. There- 
fore, an approximate solution to these equations will be 
obtained by means of inner and outer asymptotic series ex- 
Dansions for both ~ and 8. The inner expansions will approxi- 
mate local behavior near the surface of the plate in the 
Wicrmity of the discontinuity, and the outer expansions will 
model velocity and temperature profiles at greater distances 


mom the plate. 


C. BOUNDARY AND INITIAL CONDITIONS 

The nondimensional form of the governing equations given 
meow and (7) must satisfy the original physical boundary 
conditions as specified in (3a) which reduce to the following 
conditions in terms of the stream function wW and temperature 


eestribution 0: 


rye = 0: u* = sor = 0 
v* = “3 = eee ee 
vt = “3t = 0, -L<x<0 (8a) 
QO = 1 


ay > © h6Uu* = 0 
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where €* is related to the physical blowing or suction rate, 
Vy? by the following: 


Vie 5 
a = ere (64Gr) i E:3ep) 


Since this problem involves a perturbation of both the 
momentum and thermal boundary-layers which began developing 
memene Jeading edge, the initial conditions at the origin of 
maemaiscoOntinuity, 1.e., x = 0, must match these free con- 
vection profiles. For a semi-infinite vertical plate main- 
Pained at a constant temperature the following ordinary 


differential equations describe unperturbed free convection 


behavior: 
Fit’ + 3RR" - 2F'2 + H = Q (9a) 
+H" + 3FH' = 0 (9b) 
Pr 


with the following boundary conditions: 


Poeej= O: F(0) = F'(0) ; H(0) = 1 


AS mo: F'(o) = 0; H(~) = 0 (10) 


where F is the similarity stream function and H is the 
dimensionless temperature distribution and primes indicate 
derivatives with respect to the free convection similarity 


variable Ng defined as follows: 


_y Gr* 
ie = — ) fede) 
f ‘Ze1) | 7 





Zo 





where Gr* is evaluated at x+L according to the free convection 
length scale whichhas its origin at the leading edge. It 
should be noted that all references or comparisons to the 
unperturbed free convection problem employ this translated 
length scale, whereas the length scale for the perturbed 
problem is based on L with origin at x = 0. At x =0 equation 
(11) reduces to Ne = y*, The solution to the equations (9) 
are well known; however, an excellent summary of solutions 
for this and other more general cases can be found in Yang 
i |. 

eemce the initial length of the plate is unperturbed, the 
velocity and temperature profiles at x = 0 will be given by 


the following: 


(4Gr)% F'(y*) aloes 


nity ) (12b) 


u 


Cc 
Therefore, the initial velocity and temperature profiles 
Originating at x = 0 must match those obtained from equations 
(12). The details of expressing F and H in a form which will 
permit this matching condition to be mathematically formulated 


are presented in Appendix A. 


III. METHOD OF SOLUTION 


Pee REDUCTION OF THE GOVERNING EQUATIONS 
The nondimensional forms of the governing equations for 


mie problem under consideration are: 





oy atv gg ae Oty, Oty 

ay* oxoy* Ox dy*2 ove me 9) 
AL a A a*6 (7) 

oy OX OX ye Pr oy*2 


ime order to apply the Goldstein series method of solution to 
these equations suitable expansions of wp and @ must be con- 


memmeted Of the following form: 


1 
v= ©*(E/0;)7 + 6,6°[E<£, (0) ] (13) 
and 9 = 1+ $,€°(e%g, (n)] (14) 


where & =o X°; n = b,y*x?; the Wes are spossiply tunetions on 


e*; and the a, b, c, d are constants to be determined. The 
following conditions must be satisfied which impose certain 
Mmestrretions on the o's and exponent constants: 
1) The u velocities from the free convection region 
(x<0) and from the perturbed stream function, i.e., 
u* = ae Must amacem at x = 0. 
Bye the temperature profiles must match at x = 0. 
oS) The highest order derivative of the coefficient func- 


tions must appear in the zeroth order expansions which 


a7 





result upon substitution of equations (13) and (14) 
Imtovequationomo) and (7). 
4) The exponent constants in the expansions referred to 
Meo) above MUSt be antegers. 
The details of applying the conditions outlined above to 
determine the unknown parameters are presented in Appendix B 
along with the boundary conditions which must be satisfied 
By the coefficient functions fh and Sy. 
Using the results of Appendix B equations (13) and (14) 


are transformed into the following: 


y= Sere? + c2feXe ] (15a) 
and @ = 1 + E[€g,] (15D) 


MS and n = (fx) 1/Sya = yase Ce), 


where & = Gx) 
The momentum and energy equations (6) and (7) are transformed 
mn to : 
joer Kea: k K pr 
eee 15 fh) [6 °£2) - [Ck + 2)ef,) [E81] 
- de*eceser] = [een] + € + e2pek%g,] (17) 
and [(k + 1)é%g,](E£f] - ((k + 2)€%£,] (6<gt 


- de *E (ES! ] 5 S-[E gt] (18) 


where the repeated index k indicates the standard Einstein 
Summation. By the process of equating coefficients of like 
powers of —€ the following ordinary differential equations 


involving the coefficient functions are obtained. 
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Momentum Equation: 
iene order: i ieee ele = 
0 0 0 0 


Ist order: Ae ost meen otf =-—detf' - 7 (19) 
1 6 1 ole Sea Ain oe 0 


Miyo Order: ioe pepe ede fleet ety = 2Fte _ 5F Ft 
2 qo Qre2 0 2 1 ae 


ee oe. oS 
: 0 

Emercy Equation: 
wth. order: oe Cit eewe me Tae = ( 

is 0 0~ 0 0~ 0 
moe, Order: oo pet 2h cre, + eee ee See S| 

; (20) 
gZnd order: BEE Pye ieee Dee F 4git, - Bele. ae dice 


peepee i SE oe 
From Appendix B the following boundary conditions must be 


Satisfied by the zeroth order momentum and energy equations: 


ee n= 0: £ C0) = £'(0) = g (0) = 0 
A ae 2 5 0 b 
->0CO * = ees ee 
S nro: = as; - 


where i a (eo ne: DL = eee) 


The solutions of the zeroth order equations are now straight- 
forward and yield the following: 
= 2 
ie = ain (21a) 


and o, = CZ2hbs 


l 
lilt RMINATION OF THE HIGHER ORDER COEFFICIENT FUNCTIONS 
Mae first and higher order coefficient functions fh: and 
Sy for the stream function and temperature distribution inner 


expansions are obtained by solving the sets of coupled fifth 


Le 





@rder ordinary differential equations (19) and (20). These 
equations were obtained by equating like powers of §€ in the 
momentum and energy expansion defined in equations (17) and 
moe ithe presence of €* in the inhomogeneous parts of the 
differential equations introduces another independent para- 
meter which significantly complicates the numerical procedure 
employed in their solution. To eliminate this problem and 
isolate the effects of e* the coefficient functions can be 
redefined in terms of universal functions which depend only 
on the Prandtl number. The revised forms of fy and g, appear 
in equations (22), presented below. 


fame + ky 
ie th6 c ON 


on = + 6% 
©") ar - ae 
= + * + * 
ie Y 5 9 . Yo) z ae (22) 
Oo — ss ae * . x 2 
ae 270 = ee = ae 


oom ce tea ale AUNT Versa he huUneE1OnS oO fan adepela— 
0 10 


miememly On the Prandtl number. Since the first and higher 


~w 


where the Y 
il 


Sader equations of (19) and (20) are linear, the principle 
Se superposition applies and the equations of (22) can be 
ostmerted into (17) and (18) and separate like power of 
e* and obtain the following equations. 


Mire tirst order momentum equation becomes: 


Maly. .] - l (23a) 
i} 10 


and M [Y_ ] age e255) 
oe 1 i 0 
where M is the linear differential operator defined by: 
i 
Pee = ect = SPO Te OEY (256) 
1 bee oO 1) oO.) (U1 


il 


50 





Mieswctirst Order Energy equation becomes: 


Eee | = 0 
1 10 
prlaee [2 | = - 4Z' , 
1 1 al 00 
where E is 
1 
ee" 2f 2’ = 22 + 52" Y 
: Pr 1 0° 1 ot 00 
The second order momentum equation becomes: 
fe. | = 2y'*° = 3y y'' - Z ; 
2 20 10 PO 290 00 
ire | = 4Y' Y’ = SY Y" - 
a 21 ree Ly) ee Pa | 
and M [Y ] ae ee ey 
2 22 | Pl 2 | | 
maiere M is 
2 
eee et 2 YT = ATTY! + A ENVY 
2 2 0 2 0 2 0 


The second order energy equation becomes: 


ae 
10 
Zee 2 


10 


Elo | = 
2 20 
E[z ] 
2 
- 47Z' : 

10 
b= 
22 
1S 
F i 
a An 
2 Ea 5 


1) aes 


and E [Z 
2 11 


where E 


ZL 
1 


ja 


Let 


eau 
0 1 0 
ase) eau 
11 
sofort tO ae 


11 


ae 
10 


10 


' 
je 


l 


3Y vas 
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oi 
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See eee te ore. eS a7 
0 2 0 2 2 
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yt 
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the linear differential operator defined by: 


the linear differential operator defined by: 


Za 
10 


themibmncar dimrerential operator defined by: 


(24a) 


(24b) 


(246) 


[2a 
250) 


(256) 


(25d) 


(26a) 


(26b) 


(26c) 


(26d) 


iiestollowing boundary conditions must be satisfied by 


eoauatrons (23)-(26): 


At n= 0: ¥ = Y 
10 
' = 
10 
and Ze = Z 
10 
AS no: We & 
we 7 


li 
~< 


By 


(27) 











> eee ae > 10s vy > : 
aie 6 Jima 22 
Z 
ae = a Oe (28) 
2 O78 
Nn 
and Z 
2o +0, 2 +0, 2 #0 
ay 21 22 


moeenmecOllow directly from the coefficient function boundary 

mena. tons given in equations (B-33), (B-35) and B-536). 
Although they were not used in the present analysis, the 

mied Order universal function equations are presented in 


Pependix C. 


fees! RUCTION OF THE OUTER SERIES EXPANSIONS 

The inner series developed in the preceding section is 
Meerut Only in approximating local behavior at the surface in 
Pitmeebediity O1 the discontinuity. In order to describe the 
temperature and velocity profiles at greater distances from 
mie plate another series 1s required which will satisfy the 
asymptotic boundary conditions and also match the inner series 
at some intermediate distance from the surface. According 
memenesmethod of Goldstein [17] the outer series expansions 
for the stream function and temperature distribution should 
be Of the following form: 

bom y (YA) + By (YA) + 5 bt) (YA) # ee EM CYA) (29) 
and 


86 (¥*) + Eo (Y#) + 5 E28 (Y*) + Hes (Y*) (30) 
0 1 : 2 Oo 3 
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Smee as €>0 the original initial conditions must be satisfied, 
the following equalities must hold; thus defining the zeroth 
mieter Coefficients: 

my) = Fly”) 
and (31) 

@ (y*) = H(y*) 
mmere F and H are the free convection similarity functions for 
Seemoeream function and temperature distribution, respectively. 
mirsecOndition 1s merely a mathematical restatement of the phy- 
Sical condition which requires that the temperature and velocity 
meeriies at the discontinuity, €= 0, match those generated by 
M@ememicce CONVECTION process in the region between the leading 
maee and the discontinuity. 

Now the asymptotic outer series froms (29) and (30) are 
Substituted into the governing equations (6) and (7). Once again 
through the process of collecting terms which are coefficients 
Smee powers of €, the following set of linear ordinary dif- 
ferential equations is obtained from the momentum equation: 

DO (52a) 

eee PH yy ye = 0 (32b) 
SUI HY HSH - FPO = boyy = 8 Hy  (52c) 
Pao from the energy equation the following set of algebraic 
equations: » 

mY 8 HO oe) 


pio +y'e — o' —y 61 = 0 feel) 
) 4 bE 4 ie! | 





Ik if 1 it ib = 

ae '6 — 6 x a as '6 ap oul Q! eee Q! oe ar 8! = ee in 3 

2 ie, 3 ve A ie we Ms 2 ie 1 2 ae 0 Pi “9 oe) 
Mae details of the solution of equations (32) and (33) can be 


found in Appendix D and the results are summarized below: 





D) Saal 
1 
© = 
eo ' (34) 
me Za ¥ 
2 D 
= _ yp)! ee 
i = (OW es od vee 
2 
6 = 0 
1 
os: 
0 =r oF 558) 
2 
D 
6 F 


we - 2y*)@! 
wa 


2 
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With these coefficient functions determined, velocity and 
temperature profiles for large y* can be generated according 


to the following expansions: 











2 ; 2 
emerges FP Cy*) + ga PCy") 1° + [3F'(y*) 
| D (36) 
a a (5 Sei) ey = ale 
2 
and 
¢ ,D 
memuicy) + ao —H'(y*)e* + e(52— - Zy*)H" (y*)e" (37) 
2 Z 


fee NUMERICAL SOLUTION OF UNIVERSAL EQUATIONS 

The universal equations (23)-(26) with boundary conditions 
(27) and (28) were solved numerically with an IBM 3035 series 
digital computer. Double precision machine arithmetic was 


Meead for all computations. A least squares correction method 
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developed by Nachsteim and Swigert [20] was used as shooting 
point routine to determine the unknown starting values for 
each set of fifth order two point boundary value equations. 

A standard fourth order Runge-Kutta initial value integration 
routine was employed to tabulate values for the free convec- 
meometunctions F, F’, F", H and H'; the zeroth order inner 
series coefficient functions EO» fi» fi g and g)3 and the 


O 


Mmeeacre Order universal functions Y , Y , Y , etc. These 
1 0 Ld 20 

tabulated values were then used to generate the inner and 

outer velocity and temperature profiles as well as the local 


Nusselt number at the surface of the plate. 


oS 





IV. RESULTS 


Pee eRESENTATION OF RESULTS 
eee clOcuLy andes lemperature Profiles 
wee lineopwoecries Proriles 
The inner series expansions for velocity and 
Memperature were computed form the tabulated values of the 
universal functions which were obtained form the numerical 
solution of the fifth order ordinary differential equations 
defined by equations (23) through (25). The velocity profiles 
were obtained for specified values of e* the blowing or suction 


Darameter according to the following equation: 
Memo pe! + (Y' + e*Y' Je + (Y' + e€*Y'. + €*7Y' JE*)] (38) 
0 inp a 2.0 Pm 22 


where the fo Y YD? StCy mane ehlincETons .Of nw. In order 
to facilitate comparison and matching of the inner and outer 
profiles the nondimensional velocity u* was plotted as a func- 
tion of y* where the relationship between y™ and n is: y*=&-n. 
Numerous sets of profiles were generated for analysis and 
comparison; however, only the profiles corresponding to the 
following parameters are presented in this paper: 
1) Blowing and suction parameter, «*, of magni- 
tudes 0.05 and +0.25 where positive values cor- 
respond to suction and negative values to blowing. 


Z)) Streamwise nondimensional parameter, &, corres- 


Monee tos, Laratios of .25, .50, and 1.0. 
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Seep lots Of u* versus v* for these parameters are presented 
me Frigures 2-13 and for ease of comparison both inner and outer 
profiles are included on the same plot. 

The temperature profiles generated by the inner 
series expansions were obtained from the following equation: 


@=1+ Efg + (2 +e*Z )E + (Z ss aed Pee, Mein 
0 1 0 Ll 2.0 21 Zoe 


2) 
where the g > 2? Tess etc are functions of n; however, as 
was the case with the velocity profiles, @ was plotted as a 
function of y* with y* = E*n. The same parameters were used 
in plotting 6 versus y* and the outer profiles were included 
for comparison. These plots appear in Figures 14-25. 

bemeOuter Series Profiles 
The outer series expansions which were constructed 
meetie preceding chapter approximated the behavior of the ve- 
locity and temperature distributions at greater distances from 


the plate. The outer velocity profiles are given by equation 


(36) which is written below for convenience: 
* ' * C , tt RA 2 
te © fy *) + SES BG, ee = 


Cee 17D 2 . . 

PO) +g GE = Wey *)i6 (36) 
where F' and F'' are the first and second derivatives of the 
mecoe convection stream function evaluated at the location of 


the discontinuity, x = 0. They were numerically calculated 


for values of the free convection similarity variable which 


is identically equal to y* at x = 0. The constant a =e Oy 
and the constants C. and dD. CEcwecompuueda Numenically trom 
Semememons (D)-16) and (D-l7) respectively. The outer velocity 
profiles corresponding to the same e«* and &€ parameters as were 
mecemiied £Oor the inner profiles are shown in Figures 2 through 
ES 

The outer expansion which defines the temperature 
meottle Valid for large values of y* was given in equation (37) 


and 1s repeated below: 


c 


D 
= 1 ' 2 1-2 _ 1 3 
ee u(y) + Ta Se lg 6 {za Eyes inva ven ie (37) 





where the constants C Ds and a are defined as before, and, 
the H and H' are the free convection temperature distribution 
and its derivative evaluated at x = 0 where y* is the free 
@onvection Similarity variable. Once again the plots of 8 
Mmersus) y* for the same parameters as for the inner profiles, 
mecepresented in Figures 14 through 25. 
meee at Iranster Characteristics at the Surface 

ie Standard method of characterising the heat transfer 
meumetcemsirtace of a heated plate 1s to calculate the local 
Meecele number at the point of interest. In this particular 
maoptem the Nusselt number at any location above the disconti- 


Miaty, x>0, is given by the following: 





| (40) 
oe oy 
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which reduces to the following form when the appropriate sub- 
stitutions are made for T and y in terms of nondimensional 


Mariables: 


Gr,* 
Nu =-(7-) [g'(0) * (Z' (0) + e#Z' (0))E 


eee) wo (Ot ses e720 (0)) E>] AD 
20 21 De 


where the universal functions Z' , ae etc. were numerically 
P70 

evaluated at n= 0. The heat transfer data is presented in 
graphical form using the following nondimensional parameter 

Nu 

* 

4 
which was plotted as a function of € for the following values 
of the suction and blowing parameter: 


ee 9.05 and +0.25. 


Mine results are presented in Figures 26-29. 


Pee DisCUSSION OF RESULTS 
1. The Matching Condition 
oe (NCeemGa=roint ef Numerical Integration. 

The end-point of numerical integration is the 
value supplied to the Nachsteim and Swigert starting value 
momecaine at which it attempts to satisfy the boundary conditions 
memes @~. ithe choice of this point has an obvious direct in- 
fluence on the starting values, (ety and g (0), anid -alks oon 
the outer series constants C COLO |e D and D defined 


bol 20 
jmaetc-16) and (C-17) and which were numerically evaluated at 


oo 





m@resvalue of n chosen as the end-point. A summary of the 
effects of varying the end-point on both the starting values 
and outer series constants is presented in Table I. For each 
Of the end-point values, Mend? listed in the table, a complete 
set of tabular profiles were generated and qualitatively com- 
pared to determine the degree of matching evidenced. The 
results of this process revealed that the smaller the value of 
SE meme Detter was the overall matching of the profiles. 
However, for fea <2.0 no further improvement in the matching 
was observed. Also, the boundary conditions at n = ~ were 
Satisfied to equal degrees of accuracy for each ee selected. 
Therefore, all of the numerical results presented were based 
on an end-point value of (het ae 2.0. It was also observed 
that the outer profiles were quite insensitive to variations 
in the constants C .? C ete. anduthatsbOr ad. £hxed. 6" sand 
—€ the parameter Nu(Gr/4) 71/4 was insensitive to changes in 
the initial values which were caused by changing eerie 
bee ine Inner-Outer Profile Match Point 

In addition to the qualitative numerical procedure 
outlined above to effect a matching of the inner and outer 
profiles, another approach was undertaken based on the appli- 
cation of an overall energy balance to the problem under con- 
Sraeration. The principle of conservation of energy requires 
that all energy entering a specified control volume must leave 
that volume providing there 1s no storage of energy within the 


volume. In the case of the present problem all of the energy 
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TABLE  I[ 


SUMMARY OF STARTING VALUES AND MATCHING CONSTANTS 


“end 0 5.10 5.0 7.0 10.0 
Y" (0)| 2.5787 2.8611 3.0893 3.1745 
Z* (0)) -1.5531 Ries | Siesze 1-1 4004 0403 
© et aa 
Be 0) 2.5437 -3.3282 | -4.0277 | -4.3080 -4.4860 
z' (0) .12905 21959 | .345070 4034 44235 
ato 2.15407 2.3655 2.7645 2.9743 3.1229 
eos -.9905 |-1.18464 -1.4087 
= 1.605 1.636 | 1.622 | 1.615 ee 
Ce ES aE) Ce 
4.403 7.022 9.875 14.289 
-2.8477 -2.9755 | -2.9811 | -2.8945 -2.788 
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pmeecting the system at the surface of the plate must exit in 
mimuoward direction parallel to the plate. The mathematical 


moomitibation of this requirement 1s presented beiow in integral 


fo FM. 
a a yee x+L “(Dar i. 
ieee? (P - Ti judy = A oI oc) Gb (42) 
a es O O 


Mitemleft hand integral must be divided into two separate in- 
Mmeemars Since there €xists an inner and outer region where 
@iemvelocity and temperature behavior are described by differ- 
emewexpansions--the inner and outer profiles. The dividing 
maumeminctween these two regions must, therefore, define the 
femeeaepolnt of the inner and outer profiles. The details of 
eva tating the integral equation (42) are presented in 

Appendix E. The only quantative result of the energy balance 
analysis was that an exact balance could be achieved only if 
the matching point between the inner and outer expansions is 
taken to be y* = 0. That is the integral in (42) must be 
—ewlated USing the outer series to account for energy leaving 
mie volume ---the left side--and the inner series to account 

for the thermal energy entering the volume at the surface of 
Mmmemelace. AS the profiles themselves suggest, the outer 
m@rortiles when viewed in an integrated or global sense more 
nearly satisfy the energy balance, although in a pointwise sense 
they clearly do not, especially at the surface of the plate. 
Serec behavior at the surface of the plate whether with respect 


meme at transfer or shear stress is of primary interest, from 


42 


an engineering standpoint, knowledge of the exact location of 
mieematch point 1S not essential, since the inner series ex- 
Mamotons are always used to furnish results at the surface. 
mee velocity Profiles 

Before analyzing the velocity profiles some method of 
Characterizing the relative magnitude of the blowing or suction 
velocity should be established. To accomplish eles is i ne 8 = 
Semoared to the maximum value of the u-component of velocity 
in the case of unperturbed free convection evaluated at the 


location of the discontinuity, x = 0, which will be designated 


Dy U. 


. The ratiov /u which characterises this comparison 
max me oniax : 


is reduced to the following nondimensional form through 


mepropriate substitutions: 


Vv 
W 


u 
maxX 





a 
Slee GE (43) 


For a Gr = 10° for which the laminar flow regime is still 


iomeamle the relationship expressed in (43) becomes: 


=i 


WwW 





= OES) a tg (44) 


G | 


Max 


From this correlation the blowing and suction parameters used 
mieene present problem correspond to the following specific 
relationships: 

e* =0.05 implies Vv = .00256 uo, 


and (45) 
e* =0.25 implies Vv FOL 28. a 
W max 


I 





These values for ve are consistent with the original assumption 
underlying the governing boundary layer equations. 

A careful review of the inner profiles plotted in 
Figures 2-13 reveals the following features about velocity be- 
havior with suction or blowing: 


1) As suction is increased at a specified value 


OL = the location of the maximum u* is drawn 


in closer to the plate and its magnitude is 
decreased. 

2) At a given value of suction the location of 
the maximum u* is moved closer to the plate 
and the magnitude of u* is decreased as the 


value of = 1s increased. 
3) The results for blowing exhibit the reverse 
ertects on the location and magnitude of u* 
asvthese deserined in 1) and 2) “for suction: 
Numerical data on the effects of suction and blowing on 
the location and magnitude of the maximum u* for the inner 
velocity profiles as well the corresponding data for the 
Unperturbed free convection case is presented in Table II. 
Comparison of the suction data with the free convection 


mata reveals that the location of the maximum u* is shifted 


closer to the plate and also that its magnitude is reduced by 


an 
A | 


15% when ~— =0.25 and by 37% when es 1.0 with a suction mag- 
nitude of ¢* =0.05. When the suction is increased to ¢*=9.25 


the effects are correspondingly amplified. Although this is 
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mae behavior which would be expected for suction, the percent 
meduction in maximum velocity predicted by the inner expansion 
meeprobably too great considering the relative magnitude of 
Suction applied. In the case of blowing the inner series 
Mmpmerse tO produce results consistent with expected physical 
behavior since neither the location of the maximum u* is moved 
farther from the plate nor is the magnitude of u* increased 
above that observed for free convection. Also the degree of 
fmatening between the inner and outer profiles lessens as the 
metion magonitude or distance from the discontinuity is 
increased. 

As indicated above the results obtained in the case 
SG SUGCtiOn are consistent with the physical situation since 
in this case heated fluid is drawn closer to the plate where 
Viscous forces tend to reduce the maximum velocity in the 
Sereamwise direction. In the case of blowing heated fluid 
Meepushed farther from the surface away from the influence 
Mmemviscous forces and into the region where buoyancy forces 
act to accelerate the flow in the streamwise direction there- 
By increasing the maximum u*. In addition for a specific 
value of either suction or blowing the effects described 
above will be accentuated as the distance from the disconti- 


Marty 1s increased; i.e.; for increasing values of Owe ig 3 


ce 
iB 
Since the blowing or suction effects will have had a greater 


opportunity to predominate over unperturbed free convection. 
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ime lack of positive correlation in the case of blowing versus 
Mieerturbed free convection might be attributed to any of the 
mollowing: 

1) Insufficient terms in the inner series expan- 
sion; since the signs of the coefficient 
functions alternate in a pairwise manner, 
addition of another term would possible serve 
to effect a better correlation. 

2) The etfects of blowing become significant in 
thesreelon where the inner Series no Longer 
provides a valid representation of the velocity 
behavior. 

The effects of suction and blowing on the outer velocity 
foros are graphically presented in Figures 2-13; however, 
Mi@emsopecitic effects on the location and magnitude of the 
maximum value of u* are given in Table III using the same 
fomitaweas table Il for the inner profiles. The results in 
fmtewcase Of suction are somewhat inconsistent with the behavior 
mag@aeecaerelative to the free convection case. The location 
meecieemaximum velocity does not move closer to the plate, 
and its magnitude does not decrease when compared to the free 
convection profiles. This anomaly suggests that the effects 
Miesuction on velocity are more significant in the region 
Eitaseeto the surface of the plate where the behavior is better 
approximated by the inner series expansion. Conversely, in 


the blowing case the results correlate well with expected 
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TABLE III 


moCA! TONSOE SOUTER SEROPILE MAXIMUM VELOCITIES 


noi 


Nc 
TL: 


ae 





Free Convection: 
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memavior. Relative to the free convection case the maximum 
mole Of u* moves away from the plate, and its magnitude in- 
@reeases from 1.0% to 10.8%, with a blowing parameter e*=-.05, 
as x increases from .25 to 1.0. When the blowlag magnitude 
mee = -0.25, the value of een IMeGreaseset rom 2,75 todo. 2% 


as immereases from .25 to 1.0. Whether or not this increase 


ued ba (es oe 


/memex cessive must await comparison with experimentally obtained 
data. In any case it appears that a very small perturbation 
of the momentum boundary layer at the wall, produces relatively 
Peomeeticant changes at increased distances from the plate. 
ee lemperature Profiles 

The effects of blowing and suction on the temperature 
faemebes are not as visible in the plots Figures 14-25, as 
mas tme case for velocity since there are no maximum velocity 
peaks to observe. However, the same general pattern emerges 
as did with the velocity profiles, and key points are summarized 
mewovweror the inner temperature profiles: 

iim 'Ou sa given locatton above the discontinuity the 
PUNeCmeeTOLI esol VOLS, ClOoGkWwESe a2aD0UG @ =e 1.) as 
tiemMaGiitUde. Of SUGEION! 15, 1ncreased, this 1e- 
presents an increase in the gradient which is 
characteristic of an increased heat transfer rate. 

2) For given magnitude of suction there is a marked 
increase in gradient as the distance from the 
drPScOnEIMuUney In Ene streamwise direction 21S in- 


é », < , 
Greascd a,c. tne Value, or yr increases. 
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Figure 14: Suction Temperature Profiles: e¢* = 0.05, 
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Figure 15: Suction Temperature Profiles: e* = 0.05, 
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Sliettonslempera ture ;Protiles.: 
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Figure 17: Suction Temperature Profiles: e* = 0.25, — = 0.25 
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Peeure 18: Suction Temperature Profiles: e* = 0.25, a =:025U 
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Figure 20: Blowing Temperature Profiles: e* = -0.05, 
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5) Relative to the unperturbed free convection profiles 


x 


1 
7 1S to 


the effect of suction for a given value of 
increase the gradient. 

4) For blowing the reverse effects of those outlined 
in 1) through 3) above are observed with the addi- 
tional feature that the inner profiles are bowed 
Cemacuwelt deereasoumy tiieir Oradient, This be- 
havlomeehakaetenizZes sa deereased heat transter 
rate. 

These observed phenomenon correspond with what physical 
mmtuition suggests should occur. In the case of suction heated 
fluid 1s moved closer to the wall thereby reducing the thick- 
ness of the thermal boundary-layer and enhancing the rate of 
heat transfer aS 1S evidenced by the increase in gradient. The 
SepOsite Situation in which blowing forces heated fluid 
farther from the plate causes the thermal boundary-layer to 
thicken. This creates an insulating effect which serves to 
reduce the rate of heat transfer evidenced by a decrease in 
gradient. 

The outer temperature profiles again reflect the 
Mmmwemictive Shape of unperturbed free convection profiles. 

Im the case of suction their behavior is inconsistent and 
mieymdo not behave relative to the free convection profiles 
Pomde the inner profiles. This again suggests that the effects 
Of Suction are most evident near the surface of the plate 


where the inner series is a closer approximation of the actual 
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behavior. For the case of blowing the outer profiles consis- 
Bently move away from the plate with increases in either 
Maemictude Or distance above the discontinuity. Also, rela- 
Mmmemeo the free convection profiles, they behave as expected 
Mleving tarther out from the surface as the magnitude of 
Dlowing is increased for a specified value of a This suggests 
mate the outer profiles are indeed a good approximation of 
Memperacure behavior at large distances from the plate. 
Wee neat Transter at the Surtace of the Plate 
inwerdemeuomamalyze themeffeeces Of discontinuous suc- 
tion and blowing on the rate of heat transfer at the plate a 
graphical comparison was made of the following three physical 
Situations: 

1) Unperturbed free convection from an isothermal 
prbate . 

2) Free convection from an isothermal plate with 
continuous suction or blowing applied over the 
entire length. (from the study by Sparrow and 
Gess oly. 

Weeeree CONMmec don from an isothermal plate with 
discontinuous suction or blowing applied at a 
dicswamec= apove ene leading edge: 

For each of these cases the heat transfer data was correlated 
Mecording to following nondimensional parameter: 


Nu (Gr/4) % (46) 
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where the Nusselt number in each case Nu , Nu, and Nu was 
0 l 
evaluated according to length scale associated with the 


imeseontinuity; 1.€., L. The three parameters have the follow- 


me LOrm: 
1) Nu (Gr/4)~* 5 2 jf" (0) Coenen e (47) 
2} Nu, (Gr/4) Ss] HY(O)(x + 1)7% - 1.40e% (48) 
3) Nu(Gr/4)"* = - [H'(0) + (Z' (0) + AZ" (O)E 
+ (2" (9) + e*Z" (0) st e*°2' (0))E*] (49) 


The heat transfer parameters were plotted as functions of x, 


ieee. , - Hon SUcCeETONand Dlowing magnitudes, e€*, of .05 and 

.25 in Figures (26)-(29). As expected in the region from the 
leading edge, x = - L, to the discontinuity at x = 0 the heat 
transfer rates for cases 1) and 3) are identical and in case 

2) the rate is higher by a fixed amount over the unperturbed 
free convection case. At the discontinuity the curve for the 
discontinuous case experiences an abrupt change in slope and 


: : ‘ De 
intersects the continuous case curve at approximately ean 


0.22 for «* = 0.05 and at = 0.52 for e* = 0.25. The magni- 


tude of applied suction apparently has a significant effect 


on the location of the crossover point where the discontinuous 
heat transfer rate exceeds that for the continuous case. The 


behavior of the discontinuous heat transfer rate as ~ increases 
does not appear to approach an asymptotic suction limit charac- 
teristic of a constant boundary-layer thickness. However, 


Since the inner temperature Series expansion constructed to 


approximate behavior near the surface is only valid in the 
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Figure 26: Comparison of Heat Transfer at the Surface of the 
BPleateswirth Suctaon, 6° = 0.05, 
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Figure 27: Comparison of Heat Transfer at the Surface of the 
PEateownal SUucCEMon. c= 0,25, 
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Figure 28: Comparison of Heat Transfer at the Surface of the 
Plate with Blowing, €* = -0.05. 
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Comparison of Heat Transfer at the Surface of the 
Plate with Blowing, e«* =-0,25, 


80 





Oo a) 





Meomehiiborhnood Of the discontinuity, i1.e., small values of rs 


it cannot be expected to accurately predict results in asymp- 
merac limiting case as becomes large. The addition of higher 
terms to the inner expansion would serve to extend the range 

of validity, however, since the Sparrow and Cess analysis of 
the continuous case employed only first order terms, it is 
reasonable to expect the results of the present analysis to 
Meovide a closer approximation to the physical situation. 

This essential difference between the two studies should be 
kept in mind when comparing results. The situation with 
Pihowineg applied as depicted in Figures (28) and (29) shows the 
expected reduction in heat transfer rates when compared to 

the unperturbed free convection in both the continuous and 
discontinuous cases. At the discontinuity, however, the curve 
representing discontinuous blowing exhibits a downward shift 

in slope and then follows a path which is basically parallel 

to the free convection and continuous blowing curves. The 
Sm@ssover that occurs in the case of a blowing magnitude of 

e* = -0.05, is probably an anomaly due to an insufficient 
number of terms in both the continuous and discontinuous series 
eepeansions. For a blowing parameter of e* = -0.25 all three 
curves appear to asymptotically approach the same slope. A 
final judgement as to the degree of accuracy with which the 
present analysis models the heat transfer behavior with discon- 
tinuous blowingor suction applied must await comparison with 


experimental data. 
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Vi CONG@mUS TONS AND RECOMMENDAT LONS 


An overall review of the results obtained from this analysis 
Meads tO a number of general conclusions. An apparent range of 
validity for the expansions can be estimated from a comparison 
of the inner and outer profiles presented in Figures 2-25 based 
on the degree of correlation between them. In terms of the 
streamwise distance from the discontinuity a qualitative upper 
bound of : = 1.0 1s suggested. In terms of the blowing and 
suction parameter the relationship given on equation (44) sug- 
gests an upper limit of |e*|= 1.0 which equates to a magnitude 
Sieapproximately .05 Une Ons vi, homevnlch stiesOricinal boundary, 
Mayer equations should remain valid. In the case of suction 
the inner profiles appeared to model the expected behavior of 
the location and magnitude of the maximum velocity with respect 
to the unperturbed free convection pvrofiles better than did 
the outer profiles. When blowing was applied, the outer pro- 
files correlated with anticipated behavior in this respcet 
better than the inner profiles. In general the matching of 
the inner and outer velocity profiles improves as the stream- 
wise distance from the discontinuity decreases and is also 
qualitatively better in the blowing cases. The analysis also 
indicated that while the inner series provides reasonably good 
results at very small distances from the plate, in an overall 


integrated sense the outer profiles exhibit a closer correlation 
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mith the velocity behavior associated with free convection. 
The effects of suction and blowing on the temperature 
fistribution both appear to be well modeled by the inner 
serles expansions. Again the degree of matching between the 
inner and outer profiles significantly improves as the stream- 
wise distance from the discontinuity decreases. Also the 
effects of the discontinuity on temperature behavior, relative 
to the free convection case, are less direct than on velocity 
behavior. This difference is evidenced in the manner of coup- 
ling exhibited by the ordinary differential equations which 
define the inner expansion coefficient functions (see equations 
meomeand (2Z0)). This fact in addition to the apparent over- 
estimate of velocity magnitudes produced by the outer Series; 
the underestimate by the inner series; and the heat transfer 
behavior predicted for small suction magnitudes suggest that 
additional higher order terms be added to both inner and outer 
expansions to achieve a better approximation of the physical 
Situation. In general the results obtained for the effects 
of discontinuous suction or blowing on heat transfer at the 
aiiirrmace Ot the plate seemed consistent with the results for 
unperturbed free convection and the case of continuous suction 
or blowing applied over the entire plate length. Comparison 
with experimental data when it becomes available will enable 
a more definitive conclusion to be made regarding the validity 


Mameois SOlUutiOnN with respect to heat transier performance. 
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Also, with the addition of higher order terms to the 
inner series it would be possible to determine if the heat 


transfer results predicted by the present analysis tend to 
approach the asymptotic limits for large values of . SS 


tablished by Merkin [9]. In that study the following 
es mpcotic limits for continuous suction and blowing at 
mieee distance from the leading edge were reported, given 


below in terms of the parameters used in this analysis: 


= Gis PA fe xX 

1) For suction Nu/ (77) encom. = 
; Gr,% % 

2) For blowing Nu/(7-)" > 0 as Tee 


The asymptotic limit in the case of suction 1s represented 
by the dashed line in Figures (26) and (27). 

Although an exact determination of the inner-outer series 
match point is not essential when the primary interest of 
the investigation is to predict behavior at the surface of 
the plate, a more rigorous determination of that point would 
serve to enhance the overall creditability of the present 
method of solution. Therefore, another recommendation to 
improve the method point is to use the energy balance 
integral given in (42) to form an energy difference, i.e., 
from a difference integral by substracting the right side 
from the left side. <A numerical minimization procedure 
could then be applied to this difference to determine the 


match point which most closely satisfies the energy balance. 
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APPENDIX A 


COVSTRUCUPON TOP the FREE CONVECTION MATCHING SERIES 

The equations which define the behavior of the free con- 
vection velocity and temperature profiles were presented in 
equations (9a) and (9b), with appropriate boundary conditions, 
iieeerms Of Similarity functions F and H which characterize 
the stream function and temperature distribution respectively. 
In the region near the surface of the plate, y* = 0, and at 
mes location of the discontinuity, x = 0, these functions can 
Seman proximated by Suitable series expansions in powers of 
veowas tollows: 


Bey") = 2 


J 


a, (y*)? CAcA) 
O 


H(y*) = 1+ 5 b.(y*)? (A- 2) 
eee 
mieme the variable y* can be used Since it 1s identically 
equa to the free convection similarity variable at x = 0. 
The constant coefficients in (A-1l) and (A-2) can be 
determined from applying the following previously enumerated 
boundary conditions: 
At y* = 0: F(0) = 0 
Fee) = 0 
H(0) = 1 
(A= 3) 
aoe yo > 3: Fi() = 0 


H(~) = 0 
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These conditions immediately yield the result that a =a = 
0 1 
By evaluating the series representations for F" and H' the 


constants a. and DL are specified as: 
a= eae (105), ijl DT = H'(0) 


where both of these initial values are well known and depend 
@mly the Prandtl number. 


Memevaluate the higher order coefficients the form for 


mee emust be obtained from equation (9a): 
2 
Eee = = SRE’ + 2F' = A 
marem wien evaluated at y* = 0 yields F'''(0) = - 1 and 
equating this value to the series form of F’'' yields: 
m. - i 
3 6 


A similar technique is used in determining b from equation 
2 


(9b) and the following value is obtained: 


b = 0 
2 
The remainder of the coefficients are found by successively 
Gitrerentiating equations (9a) and (9b) and equating them to 


the appropriate series representation evaluated at y* = 0. 


The results of this process are summarized below: 


a = 0 b = H'(0) 
0 1 
a = 0 b = 0 
1 2 
(A-4) 
a, = 2 (0) D, = 0 


9 
i} 


eVebre (0 He 0) 
- 1/6 b = we oe 
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-H'! (0) 


a — 
‘ 24 
tt 2 
a Aaa) eh 
2 120 


itese series forms for F and H will be used to satisfy the 
matching conditions for both the inner and outer expansions 


later in the course of the analysis. 
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Pelee eNO Y B 


DETERMINATION OF THE STREAM FUNCTION AND TEMPERATURE DISTRI - 
BUTION SERIES EXPANSION PARAMETERS 


Recalling that the definitions of the stream function and 
temperature profile series expansions involved certain arbi- 
trary constants, the methods and procedures used in evaluating 


eeemeare Outlined in this appendix. 


l 


b= e8(E/o J? + 6 ET [ESE Cn) ] (B-1) 

a= 1+ 6 6° re%g, (n)] (B-2) 
where 

E = x* and n = by ax? (B-3) 


The object of constructing such expansions for the stream 
Miterton and temperature distribution is to obtain series re- 
presentations which satisfy the momentum and energy equations 
which mathematically model the physical problem. The indivi- 
dual components of those governing equations (6) and (7), l.e., 
the partial derivatives, must be evaluated in terms of 
Semesstons (B-1) through (B-3). Substituting (B-1) and (B-2) 
into the momentum and energy equations yields: 


Momentum Equation 


1-2b 2b-1 
so 9 920 6° * 8 ffate + k) + been] (e*ey] 
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| bo 
ey 
i) 
oO 


Emerey Equation 


isle Apel 
a a ee k key 
49 oo 6 8 {[a(d + k)Eg [6° £71 
k k 2 2+ d 
nice mom e se Io 4a" 6 E> Ey 
ts Z 
Loe, td ye 
= pe o'o [egy] (B-5) 


Certain conditions which must be satisfied will impose 
restrictions on the allowable values for the o's and constants 
a, b, c and d. The one condition which must be satisfied by 
the equations 1s the initial velocity and temperature matching 
at a distance L from the leading edge of the plate where the 
discontinuity in normal velocity begins. At this location the 
perturbed velocity and temperature profiles must match the 
free convection profiles. In terms of the stream function the 
memecity 1S given by: 

ut = Bee eho Po SP poke 27 EP ey B68) 
The free convection velocity profile written in terms of u* 


GCanmebe expressed as follows: 


u* = F'{v*] = J yee Gos (B-7) 
aD 


J 
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Pawating the leading term in each series (B-6) and (B-7) i.e., 


when j=2 and k=0 we obtain the following relationship: 
fed aco | 

Pi 1 

2ay* =o > g ZH? Ff, (B-8) 

2 12 oe 


For proper matching, the exponents of y* must be equal; therefore, 


ac z : 
po - 1-1 (329) 


Puescr, Since it is desirable to have expressions which are in- 
dependent of ¢* for the boundary conditions on the fy as nr, 
the following condition will be imposed: 


ca 


oo > 6 = 1 (B-10) 


1 2 3 
The requirements formatching the temperature profiles yields: 
_ad _ka _ka 
fanee, P k 
o (—) [oro 
4 TT ie 


. | 
P(E)? gpd = jE bay”? (B-11) 


Applying similar arguments to those used for the condition of 


aoe. 
1 2 


velocity matching, the following relationships are obtained 


meom (B-11): 
ad _ - 
i il (B-12) 
and 
696 6 = 1 (B-13) 
oe 2 4 


femmeacgditional restrictions on the exponent constants consider 
equation (B-4) in which the highest order derivative term must 
appear in the Oth order expansion in order to obtaine meaningful 
results. This requirement imposes the following condition on 


the exponents of &: 


c + Sb. 2c + 2b-1 


a ae a 
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Sreecaquivalently, 


ae = [5 ea A Ik (B-14) 


Mitemsame restriction also applies to the energy expansion, 
fauation {B-5) and the same relationship is obtained (B-14). 
Since fractional exponents on § are not permissable, due to 
meememability to match coefficients of like of powers of € in 
both expansions, from equation (B-12) the following further 


[irormation is obtained: 
d= 1 [B= kS:) 
and 


= ] (B-16) 


oo 


Now combining (B-9) and (B-16), the constant cis specified 
as: 


c= 2 (B-17) 


UWsame this information and the relationship given in (B-14) 
imposes the following restrictions on the possible values of 
a and b: 

meal) OS and b < = 1/3 (B-18) 
iiemectlationships previously obtained for the o's can be sim- 
memered bY Substituting the known values for the exponents in 


equations (B-10) and (B-13) from which are obtained the 


following: 

oro > = (B-19) 
and 

6-99, = 1 (B-20) 


eal 





Samece Cquation (8-5) is homogeneous in om and on appears only 

Smelly in equation (B-4), without loss of generality it can 

be set equal to unity. Setting >. = 1] in equations (B-19) and 
(b=20) yields the following further relationships for the 6's: 


oe = il (B-21) 
1 2 
and 
oy = 1] (Ba22,) 


The foregoing analysis has permitted the reformulation of 
equations (B-4) and (B-5) in terms of only two unspecified 
parameters a and > > and those equations reduce to the follow- 
mc £Orms : 1 1 


— = — 


Momentum: dag °E frtk + L) eset] (eS€4] Sle at ae es 
[eserj} - de*eeXey] = (eSeue] + € +e7[e%g,]  (B-23) 


= 
: . ary k ke, k 
nergy: fad E Clee lise g  ILé fy] Sil GK ees) Gell 


m |e 


(Eat ]} - 4e*E(ESgt] = selé et] (B-24) 
Equation (B-18) allows for only two reasonable choices for 
meeevailue of a; namely, a = 1/3 or a = 1/2. Furthermore, for 
each of these two values >, can be defined to either include 
= Or not; thereby posing four possible cases for consideration. 
The cases which specify a = 1/2 do not satisfy the boundary 
Senaition for fo as n>: and therefore, are unacceptable for 
memo Gesent problem. The case for which ¢* is included in ve 


and a = 1/3 results in a set of ordinary differential equations 


a2 





which include e«* in the homogeneous part and which cannot be 
transformed into general equations independent of e¢* through 
the use of universal functions. Therefore, the case with 
mel, > and €* not included in the definition of > was 
PeereGted tO complete the solution of the problem. A summary 


of the results is presented below: 


meee /S. b = = 1/3; c=2, d= 1 (B-25) 
3 IGS Be WS 
oe zi (Z) 9 o. ~ (7) ’ oe = 2 =k (B= 204) 
es ob 5S 
5 5 * 
E= (7) ,n= Ge y*=%/, (B-27) 
b= pete? + er feXe | (B-28) 
k 
me eee | (Ba Z9) 


The momentum and energy expansions which result from these 
definitions were presented in Chapter II, equations (17) 
and (18). All that remains is to specify the appropriate 


boundary conditions for the fy and Sy: 


5 k . 
At n=0: u* = ae = ELE%£t] = 0, (B-30) 
ani 
ve = = fee + Le [cKeayeXe,]} = te*, (B31) 
and 
@= 1+ Efe%g,] = 1 (32) 


From these original boundary conditions, the local boundary 
conditions for fy and fy are determined to be: 
= ' = = _22 
£, (0) £,. (0) g, (0) Oiee bows sa lek (B-33) 


The boundary conditions at n-~ are obtained through matching 


@iesvelOcity and temperature profiles to the free convection 


oe 





profiles. The velocity match is performed by equating equations 
feo) and (B-7) with the appropriate values substituted for the 
unspecified parameters. The result of this process is the 


moneowing relationship: 


Rt) a eer on = 
j22 ja, *) [ee Eg (B-34) 


By equating like powers of y* and substituting the values for 


INE 2 obtained in Appendix A the following local boundary con- 


Gace aons for the fy as n>o are determined: 





_' 
— = 2a = F"(0) as nr 

1 2 
nN 
a 
i = Se) = 2 Je as n7o 
Ae 3 2 

(B-35) 

' 
4 =da =- aD, as n>@ 
n3 ly 
ia ‘< 2 
—2 = Sq = IES Cee (0) J as nro 
n* 5 24 


iMiemecemperature profile matching yields the following local 


boundary conditions for the ere 


S 

_/ = b = eel OR as nwo 

S 

a2 = b. =a) as nro 

. (Bi=.560)) 
aS = b. =i as nro 

g , PrF"(0)H' (0) 

ne 2 = oe as n7o 
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UE eNO AG 


THIRD ORDER UNIVERSAL FUNCTION EQUATIONS 
The third order ordinary differential equations which 
mMeoult trom the expansions given in (1/7) and {18) are pre- 
sented below, although only functions of second order 


and lower were employed in the actual analysis: 


Peewee fet OL fee Sf E' = SEI f= giver 
3 0 3 0 3 0 3 a ey) t22 iene? 


Ee de xf" - g (ee 
2 i 
1 
—g" + 2£ g' - 4£'g + Sg'f - g £' = 2g £' + 3g Ff 
Pres : Vee See 5 3 E 3 eee ee 
= t_ t 2 ae - 
3f 8. ee Neo (G=2) 


itwemGorresponding forms for f and g in terms of the univer- 
3 | 


Beene t ions and the resulting differential equations are as 


follows: 

f = ee + a a + ray, + a (C232) 
and 

g= Z ee ert ect a ete eZ (C-3b) 


3 3 0 oY Che es. 


The third order momentum equation becomes: 


Gummo Ya OY YM. = Ay" Y SA ie (C-4a) 
3° 30 10 20 Vices aig 150 2210 1 0 


Mm [Y..] Ser emer ear SY yi te SY 8 YN 
| ees Peo ioe u f0n a 11 20 


o5 





nee = 2. (C-4b) 


LiOae 25 2 0 a 
M (Y J= S¥y' Y' + 5Y' Yt - 3Y Y™ - 3y y" 
3 a2 0 22.2 Es se? a | eG: 2:2 as ee | 
eAyUNY) =04y' Y = 4". (C-4c) 
105 2-2 i ee a 2) 
Py je= SY’ ¥' = SY YY" = ayy = ay (C- 4d) 
3 ore i ee aad 1 i aaa oe sR ea a 22 


Pee tne linear differential operator M is defined by: 
3 


eee” Se gen = Ot at st Sy (C-4e) 
3 0 3 Nae 0 


3 3 


The third order energy equation becomes: 


Emi | = 22 6% + 32° Y' = 3Y 2%) --4¥ 7°; (C-Sa) 
| 3° 1¢@ 26 220 Loo 10-206 2505 6 
Pe] 22 Ve 22s P32" a 3Z. OY 
3 Soa VO? 20 Pele Bez 20 2c zo 0 
BC EC OSV dV 7k SA 1 Ca eG Gen) 
Poe 2-31 ely 20 ZO 7~sls1 2 1 10 20 


E(Z )=2Z Y' + 22 
1st aay ae 


Ne Ee Se ee 
oar 3 2 a 


pA Zol- Alek Ze 0 


= SY ess A) CO San Mg! 9) CS ee Me) Ge 


t 
G22 Pte ee Al Zk tie a ee 


RTA ae Ae 5G.) 
Za 


moe. | = 22 Y' + 32 YO PS wee, ee ee 
1 1 2 


3 33 Sy | oa 22 Paz 22 A 22 05d) 


Winere tne linear differential operator E is defined by: 
2 


ae 52" me te 7st yg (C255) 
YT 3 0 3 0 3 00 3 


3 3 00 


The boundary conditions for these equations are: 


ot (| Ue | = = Y =e i =e ier 
30 31 32 33 
Vewe= ey! = YY" = Yr" = \0-. (C-6) 
30 31 a2 33 
and 
Z =e =e Z = Q; 
3 0 31 725 33 
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Ge. 


aiieedo NOs 4 fF" (0) ]* ee 0) a ee Oy a 0 
at 24 ea 7 "32 "33 
(C-7) 
L 1 t 
and — > als OU GU, fi eee Oe 2 0 


i 8 


Meee tollow from the coefficient function boundarv conditions 


ce 32 33 


Peeeet 17, Equations (B-33), (B-35) and (B-36). 
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Debate tON OF Tb OUTER SERIES COEFFIGLENTS 

fine mpeoOcesc On deterimmingmene Coefficient functions of 
the outer series expansions defined in equations (29) and 
(30) involves the solution of the ordinary linear differential 
emmrttens (52); the evaluation of the resulting constants of 
integration; and the solution of the algebraic equations (33). 
itemerocedure 1S greatly Simplified if the outer expansions 
for wy and 8 are first transformed into series with coefficients 
Papecliueare powers series in y*, the coefficients of which are 
derived from the asymptotic forms of the fy and gi. Thais Geeh- 
nique will facilitate the evaluation of the constants of in- 
Mereration Significantly reduce the complexity of equations 
(eo) . 

The asymptotic forms of the inner series coefficients can 
be written as polynomial in n with as yet undetermined 


coefficients as follows: 


aa 2) 1" (D-1la) 
zZ 2 

m= A n° SPB nn? + Cn + D (D-1b) 
1 l 1 1 il 

fF =An* +Bn?>+Cn*+DnietE§ are} 
Z Zz Z Z Zz 

a. = by (DS Ba) 
0 if 

g = B*n + C* (D-le) 
1 1 1 

a= eon Con + D (D-1f) 
2 Ps 2 2 
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which must satisfy the asymptotic boundary conditions specified 
mee b- 55) and) (B-56). The constants are evaluated by substi- 
meme the equations (D-1) into the original differential 
Samaeions (19) and (20) and then solving the algebraic equa- 
tions which appear as coefficients of like power of n. The 


results of this process are summarized below: 





ees eon <a 
an 7 a 5 Be 0 
Dac 
= -_ * = iomeel = 
B= 0 B= 0 % za on 
e 2a x = 
= TBD C. ae Be ¢) 
2 
D = Ze* D = TBD ce = ~biCi 
l Z Ps 8a? 
2 
E = OF D* = Fe 
? 4a 2 6a 
Z 
1 


where a Say eon) a = He) alin Cc and D to be determined 
maotmapplication of the matching condition. This process 
will be outlined at the end of this appendix. 

Now the equations (D-1) with the coefficients as defined 
mime.) Can be substituted into stream function (15a) and 
memperature distribution (15b). Then if the definition of 
fee 7 © «LS invoked and all coefficients of like powers of & 


grouped together, the following expansions are obtained which 


mem alid as y* becomes large: 


b= [ayt?- yt? - thy**] + [O]E 
G 
rij ae ee et 4 2 * 3 _2 
# [C y* - ghy**]E? + (D y*]é (D-3) 


a9 














and @ = [1+ by*] + [OJE +|5-— - —-y*!|e? 
: 2 Sa° 
b D 
e eee 2 (D-4) 
6a | 


With these equations in hand the final determination of the 
meer Series coefficients can now proceed. 
Miemarrst ot the ditterential equations defined in (32a) 
memeeoedrithmic and its solution has the form: 
p = Kyp' (D-3) 
1 1 0 
Since v is the coefficient of the linear term in the outer 
Semamsion and that coefficient in (D-3) is zero, KT must also 
equal zero, and y has Ene. final form: 
yp = 0 iso) 
Using this equation in (D-2b) reduces it to logarithmic form 
heen solution for ue as before: 
y =K yp! coe 
2 fa 0 
where K ro evaliaced by equating {)-7) with the appropriate 
memenmeot the £° coefficient in (D-3). The final result for 
vw then becomes: 
S C 
wy = ] yp! (D-8) 
2 24 0 
2 
Ime last of the differential equations (32c) can be expressed 


wimiesampler form using (D-6); however, its solution requires 


a different technique. The reduced form for (32c) is: 


= oe > it pi = © + yi (D-9) 
ey 3 eee a6 0 0 
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Using the free convection stream function equation (9a) p''' 
0 
can be written as: 


ye oe sh oy we apt eee (D=10) 


0 
and by suitable algebraic manipulations in combining (D-9) 


and (D-10) the following differential form is obtained: 


ay. af’. \_, 
a Ona Die! 
(34) + 6(5 on 
0 0 
Emmarehnt forward integration yields the following general 


Bemme1on for w ; 
S| 


y = OW - 2y*p' + Kp! (D-12) 
0 0 3 0 


3 
where K is obtained as before with the aid of equation (D-3). 


Thus the particular solution Ve has the following form: 


ies OW yey + Za (D-13) 
This completes the determination of the outer expansions 
for » and @ with the exception of the constants C and ay 
Since these constants appeared originally in the asmyptotic 
forms of : and a and are determined from the matching of the 


inner and outer series, the results of equating the universal 


Mametron +Orms to the asmyptotic forms will yield expressions 


fomee, and D . The following equations result from this 
l zZ 
process: 
6 1 3 4 k fe 
mom ts == Wi + Cn - s+ €"™ = Y et (D-14) 
1 6 9 1 0 9 10 LA 
2 
eigeerOr £ : 5a y* aa n> +Dn+t cue 
2 Die 0 4a 0 20 da 
2 2 
7 sspears + ¢*27Y (Da15) 
2,0 20 ee 


Wow 





wolvine for © trom (D-14) yields the following: 
1 





Ge = € te ; (D-16) 
i Ip he, | ie! 
where 
Y iy 2s 
aa ( n oe J; 
and 
m4 
4 
= ee sue 
om ( n i aL 


PeomilLar procedure yields the following form for D obtained 














tmome ( D-15): 
fee= D + €*D +€*7D (oat 7) 
2 Z 0 ZA 1 2°2 
where 
0 e. 3 oe Ct 
= 2 ws een) - 9 
ee om 2h * te 1 ae? 
Z 2 
u 1 : 1 & Oe 
Ze 2 1 0 
ies, ta: tan 2° 
Z 2: 
and 
yG Cs 
= Zee 11 
oer ( n 4a n? 


where the value n was taken as ea the chosen end point of 
numerical integration. This completes the derivation of the 
outer series expansions, except to note that the constants 

c Poe. anid DS are evluated numerically from the 


9 ? 
10 Ii 2 0 fas 


maombated universal functions. 
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APPENDS X E 


INTEGRAL ENERGY BALANCE ANALYSIS 
In only heat transfer problem an overall thermodynamic 
emergy balance must be satisfied. For the case of the present 


problem the integral energy balance takes the following form: 


clo (Tl - Tm) udy = (**'o" (x)ax (Ee 
0 Pp 2 0 0 
The left hand side reduces to the following form with non- 
dimensional integrand: 
7 [ae 1/4 
ono, - Ce, 


fvu*edy* (E-2) 
Since there are two distinct regions over which this integral 
is defined, it must be split into two integrals separated by 
the point, Vedge where the inner and outer expansions should 
match. Substituting the appropriate forms for u* and @ as 
defined by the inner and outer series the integral (E-2) can 
be written as the sum of the following two integrals: 

Inner region y* = 0 to y* =y, 

wv" Megrnye + [et(m)g (n) + £t(n) 18? 

+ [£'(n)g (n) + g (n)£t(n) + £t(n)]e*} dy* —— (E-3a) 


— = Gr 
we vo (T,,-T..) (q— ) 


wmnere ¥ 1S the matching point, 
* 
an = = ada) = Gdn which transforms (H-S5a) into the 


following: 


eae, 
- > Gr We tr-2 t t 3 
qieaye (ts Gi f {fe + [f"g | + fils ie 


eS 





! ! t 3 _ 
[fig, + gf! - frye }dn (E-3b) 
Miteme y depends On x. 


Popemtine Oller region, y= = y to y* = », the following integral 


aoe Les : 


-O 


_—7 _T Gr % t ] 
4c,ov(T,, Gag f f{F'H + Ta 
1 D 
| peal » aes Gi Ca a Dy AEE) Gat chy (E-4) 
2 


(F'H' ie a a esi 


Al C2 


oe Ae) 
+ [e(yh - 2y*) + 
2 

The right hand integral in (E-1) also must be divided into 
two distinct regions, the first extending from the leading 
edge to the discontinuity and the second from the discontinuity 
[memeany arbitrary location above it. The following two inte- 
Oars are obained: 

franax + requ dx (E-5) 
where gti is the heat flux in the undisturbed free convection 
meoton Of the plate and the TSeehne Medustlux abovertme discon 
mamtircy, bOth evaluated at the surface of the plate. Upon 
Making suitable substitutions to non-dimensionalize the in- 


tegrands, the following result is obtained from (D-4): 


meas ole, IS ease al ; 
-4k(T,-TA(p) (gH) +z gtet eget 
+ = g'ES] eo) 
2 


where Gr is evaluated at L and the oe are evaluated at tne 
Smcace Of the plate, n = 0. 

After completing the indicated integrations in (E-3a) and 
(E-3b) and equating coefficients of like powers of € the 


following relationships are obtained: 
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nC) moma) eee Tr Cy )Hty ) (E-7) 


C &°F' (y )H(y) = [F"(0)]*y° (E-8) 
Jen coyH' (0) (@-) + £ (%) = FE)H) - Le ug) 
S g2 aes 3Pr 
D i 1 
: Tia F' Hy) egy HG) = ae HCO); (E-3) 


Giiemrelationships can be identically satisfied only at the 
surface of the plate, n= 0. From this fact it can be concluded 
that the outer profiles are the better approximation of global 
behavior and will satisfy the energy balance in an integrated 
sense, pointwise local behavior is required near the surface 

of the plate the inner expansions should be used, although 


they Cannot satisfy the integral energy balance. 
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